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Abstract 

Analytical results for the spatial dependence of the correlation functions for all meson excitations 
in perturbative Quantum Chromodynamics, the lowest order, are calculated. The meson screening 
mass is obtained as a large distance limit of the correlation function. Our analysis leads to a better 
understanding of the excitations of Quark Gluon Plasma at sufficiently large temperatures and may 
be of relevance for future numerical calculations with fully interacting Quantum Chromodynamics. 
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I. INTRODUCTION 



In this work we study mesonic correlation functions in high temperature Quantum Chro- 
modynamics (QCD) in the lowest order of the perturbative expansion. We follow the method 
developed by [l| where the simplest pseudoscalar case was calculated. Calculations are per- 
formed analytically with a special attention to the divergent parts, which need to be regu- 
larized and we adopt a Pauli-Villars ^ regularization scheme. The range of the correlation 
function is determined by the screening mass, which can be defined as the inverse screen- 
ing length characterizing the exponential fall-off of the mesonic spatial correlator. For zero 
quark mass at the high temperature T the meson and baryon screening masses approach 
their ideal gas value, 2irT and 3irT 0, 4] respectively. 

The pseudoscalar meson screening mass was calculated jl| in the noninteracting case 
m scr = 2yJn 2 T 2 + m 2 , for finite quark mass m. Calculations in the interacting Quark Gluon 
Plasma (QGP) were done in nonrelativistic QCD fl-Q for all meson channels and in the 
Hard Thermal Loop Approximation (HTL) of the QCD [8] for the pseudo scalar one. The 
last, numerical approach, requires the analytical form of the mesonic spatial correlation 
functions in the free case in order to calculate the divergent integrals. 

The main scope of this work is to provide the analytical expressions for the correlation 
functions for all meson channels to allow for future calculations of the screening masses in 
the interacting QGP. 

The importance of such evaluation in order to identify the relevant degrees of freedom in 
hot and dense QCD was first pointed out in js-ll], where the screening masses of mesons 
and nucleons were calculated on the lattice. More recent results can be found in 
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231 ] . Though the evaluation of hadronic screening masses is a major achievement of lattice 



studies of the degrees of freedom characterizing the hot QCD (the large number of lattice 
sites available along the spatial directions allows one to study the large distance behaviour 
of the correlators), analytical approaches are not so common in the literature. 
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II. MESONIC SPATIAL CORRELATION FUNCTIONS 



The lowest order quark loop contribution to the correlation function is defined by the 
expression 



G M (q) = -2N c Tj2 



d 3 p 



Tr 



M~ 



+ 



m 



M~ 



— m 



(1) 



from which we calculate spatial dependence for static plane-like perturbation 

+00 



G M (z) 
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2^ 



dq z G M (0,q z ) e 



2\ Jq z z 
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Here M=S, PS, V, PV corresponds to the scalar, pseudoscalar, vector and pseudovector 
channel (T$ = 1, Yps = 75, Ty = 7^, Tpy = 757^), <? M = (iu m , q) is the external momentum, 
cu m = 2miiiT is a boson Matsubara frequency. Summation runs over fermion Matsubara 



frequencies of quark internal momentum p M 
the trace over the spinor indices. 



Tr 
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(iu n ,p), ui n = (2n + V)mT and Tr denotes 



(3) 
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We convert the sum over the frequencies 24j] to two integrals. The first one, independent 



of T which is called a vacuum part (divergent one) and the second, T-dependent, called the 
matter part. 
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Let's start from the vacuum part first 



+ioo 



dp f d 3 p p 2 — q 2 /4 + m 2 



— ICG 

+ioo 



G v ps{q) = SzN c 



2tt J (2n) 3 [(p~q/2) 2 -m 2 ][(p + q/2) 2 ~m 2 ] 
dp f d 3 p p 2 — q 2 /A — m 2 



-lOO 

+ioo 



G v v c (q) = mN c J 



2tt J (2n) 3 [(p-q/2) 2 -m 2 ][(p + q/2) 2 -m 2 ] 
dp f d 3 p p 2 — q 2 /A — 2m 2 



(2tt) 3 [(p - q/2) 2 - m 2 ][(p + q/2) 2 - m 
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-too 
+ioo 



G mc (a) - mN I ^ [ d3p p2 ~ g2/4 + 2m ' (5) 



-too 



After a Wick rotation (p° — >■ ip 4 and g° — )■ ^4) in the Euclidean metric = + q 2 = —q 2 ) 
we can express (jSJ) with two integrals 

h(m)=8N c [ f P * 2 ) - , (6) 

J a (m, te ) = -4iV c I (2 ^4 [(pf , _ Qe/2) 2 + m 2 ][(pf , + gs/2 )2 + m2] > (7) 
where c? 4 j?e = d 3 pdp^. Both integrals are divergent and must be regularized. Following |l| 
we use Pauli-Villars regularization and define regularized functions: 

N N 

I?=J2 a MMi) and I*{q E ) = J2 a MMi,qE), (8) 

8=0 8=0 

where ao — 1, Mo = m and higher masses are large. Coefficients a, for z > fulfill the 
following equations 

N N N 

J> = 0, J>M* = 0, ... J2 a ^ (N ^ )=0 - ( 9 ) 

8=0 i=0 i=0 

After some calculations we have 
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8=0 



J 2^) = 



2tt 2 
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The vacuum part of the correlation function can be written 

-I?-q 2 E I«(q E )-im 2 lX(q E ), 

2I« + 2qll*(q E )-Am 2 I 2 R (q E ), 
2/f + 2q 2 E I 2 R (q E ) + 12m 2 I 2 R (q E ) 



G v s°{q) 

GTsiq) 

G v v ac (q) 



(11) 



The next step is to obtain a ^-dependent correlation function. Performing the Fourier 
transforms of (fTTT) we find that in the complex q E plane the function I 2 {qE) has cuts for 
imaginary q E = ik from ±2Mj to ±00 respectively. Performing the integration (there is no 
contribution from I± integral) the q 2 E I 2 gives 
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(12) 



2M, 



i=0 



where K 2 is a modified Bessel function of a second kind. The l 2 gives 
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2A1, 







i^z 2 M 2 






2' u ' 2 / 





where G 3 '® is a Meijer G function 25|. In both summations f lT2"fT3"j) the only first term 
survives and finally 



GT(z) 
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^m 2 K 2 (2mz) + ^-m 3 M vac (m, z), 
n 2 z 2ir 2 
on S/V 

c -m 2 K 2 {2mz) - ^m 3 M vac {m, z), 



n 2 z 



where M vac {m,z) = G 3 {° 3 z 2 m 2 
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The matter part of the correlation function 



is 



G™*{q) = -AiN c / ^- 



dp f d 3 p ^ r 



M 



M 



2tt J (2tt) 



ePo/T + 1 



(15) 



where Tr is evaluated in ([3]). After change of variables p — > p' — q/2, like in the vacuum 
part, and keeping the static limit (q° = 0) we evaluate the energy integral by deforming a 
contour of integration around the poles 



p° = u±(p, q) = a/tti 2 + (p± q/2) 2 



(16) 



and obtain 
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G™ u (q) = -32N C J 



d 3 p 

d 3 p 
(2^)3 

d 3 p 
(2tt)3 
d 3 p 



pq + q 2 /2 + 2m 2 ^ pq — q 2 /2 — 2m 2 



Au + pq(e u +l T + 1) \uj_pq{e"-l T + 1) 
pq + q 2 /2 pq- q 2 /2 
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(17) 



4u + pq(e u '+/ T + 1) Auj-pq(e ul -l T + 1) 
pq + q 2 /2 + 3m 2 ^ pq — q 2 /2 — 3m 2 



(2tt) 3 [Au + pq{e^+/ T + 1) 4:U-pq{e UJ -l T + 1) 

Then we simplify by changing the variables p' — > p + q/2 in the first part of the sum and 
p' — > p — q/2 in the second part and integrate over angles. 

oo 
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where a; p = \/p 2 + m 2 , p = \p\, q = \q\ and 



A x (p,q) = 4 + - In 
p 



2p — q 



A 2 (p,q) = — In 
pq 



2p + q 
2p — q 
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Since the whole q dependence is in the A\ and A 2 functions, we only need the Fourier 
transform, see (jSj), of A\ and A 2 

+00 



A 1 (p,z) = i- j ' dq + 



2p- q 



+00 

If 2 

MijP,z) = 7T / dq — In 
2ix J pq 

—00 



2p + q 
2p — g 
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In order to perform the integration we use the identity 

In 



iqz 



(20) 



2p — q 
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(21) 



e^o2 (2p + g) 2 + e 2 ' 

The function (j2ip has two cuts on the upper half plane q = ik + 2p and q = ik — 2p. Fourier 
transform of Ai is easy to obtain 



iqz 



Ai(p. z = hm — / do < 4 + — In — ^- 

^ ; e^o2n J y \ 2p (2p + g) 2 + e 2 

—00 

sin(2p2;) 2 cos(2pz) 
pz 2 z 

The second case, A2, is more complicated and we finish with two integrals 
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so that finally we have 



A 2 (p, z) = - [T(0, -2zpz) - T(0, 2zpz)] 



(25) 



where r(0, ix) is the incomplete gamma function. In order to obtain the z-dependent corre- 
lation function in the matter case we have to insert (122123ft into (fl8l) . To perform integration 
over p we can use the identity 
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for the first term of the expansion ff26l) gives 
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For each element of the sum ( )26l) (for every Z and T) we have 
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Taking into account the relation 
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we have 



M matt (m } z) = 7rmzM vac (m } z) 
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and finally after some cancellations 
-N C T 



Gv(z) 
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7TZ 2 
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The main contribution to the summations comes from two elements / = and I — — 1. In 
the 2: — >■ 00 limit all channels give the same asymptotic 

G{z) Z ^°°> C0USt c -2V*' 2 T' 2 +m 2 z = COrLSt c -m scr z ^ 

where we identified the screening mass m scr 



m scr = 2v / 7t 2 T 2 + m 2 , (35) 
which for massless quarks has a simple form m scr = 2ttT. 

III. CONCLUSIONS 

In this work we have shown that due to the specific cancellations of Meijer G functions 
the final formulas for the spatial correlation function are compact and differ only with the 
sign or some numerical constants between different channels in the weakly interacting QGP. 

The screening mass is the quantity which governs the large-distance exponential decay 
of the correlations of mesonic current operators and is of course equal, in our system, for all 



mesons and has a value m scr = 2\/tt 2 T 2 + m 2 , where m is the free quark mass. 

We have explicitly evaluated the analytical form of the spatial, z-dependent mesonic cor- 
relation functions for all mesonic channels, which is of great importance for future numerical 
calculations of the effective masses of mesons in the fully interacting QGP, especially that 
in the literature there is a discrepancy between lattice 12|-|23| and non-lattice [5|, |6|, |8[ calcu- 
lations of the mesonic screening masses. In the lattice calculations in the high temperature 
limit, screening masses approach the non interacting value from below, whereas in non lattice 
calculations they approach from above. 



In papers [5|, |6| the analytic predictions for the screening masses, related to various 
quark- ant iquark excitations at high temperatures, are determined in the nonrelativistic 3- 
dimensional QCD effective theory with the next-to-leading-order perturbative corrections. 
In paper |8[ the numerical calculation in the HTL approximation, in the pseudo scalar 
channel only, was performed with the help of a new technique of numerical regularization 
of divergent integrals with the usage of the analytical formulas of the spatial correlation 
functions from the non interacting case, which we provide in this work. 
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